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1. Introduction 

We have witnessed a proliferate growth in theoretical and experimental efforts to 
understand and control physical systems in a quantum level. In particular, possibilities 
to apply quantum mechanics toward the radical improvement of information technology 
have magnetized attention from various branches of physics and beyond. Ever since the 
advent of lasers, optics has been closely related to the foundations and applications of 
quantum physics because of its controllability to a very fine level. However, photons 
do not interact each other, which makes it difficult to engineer their states. In order to 
overcome this problem, various schemes have been suggested. 

In a quantum optics laboratory, Gaussian states, whose phase properties are 
described by Gaussian probability-like functions, were generated but there was some 
limitation to use them for various tasks of quantum information processing. There 
have been suggestions and realisations to engineer the quantum state by subtracting or 
adding photons from/to a Gaussian field. Photon subtraction and addition to a given 
state are plausible ways to manipulate a quantum state. 

Many of the quantum optics textbooks start with the definition of annihilation a 
and creation a) operators. When they are applied to \n) which represents a state with 
n number of photons: 

a\n) = \fn\n — 1) and aJ\n) = \/n + \\n + 1). (1) 

It is easy to recognize that by annihilating (creating) a photon the photon number 
changes to n — 1 in + 1) but it is not straightforward to acknowledge the coefficient \fn 

(v^TI). 

In 1924, Bose pQ found that the most chaotic light bears the spectrum of Planck's 
formula at a given temperature under the two assumptions 1) light is composed of 
indistinguishable particles and 2) any number of particles can occupy one quantum 
state. This is the first full quantum-mechanical explanation of the Planck's formula for 
the blackbody radiation. Einstein extended Bose's idea to the quantum theory of an 
ideal gas pj. As a whole, the particles which bear the statistics under Bose's assumptions 
are now called as 'bosons'. Bose's first assumption is against his contemporary belief 
that particles are always distinguishable. In the theory of identical particles, Bose's 
second assumption naturally brings about the symmetric nature for the state of bosons 
in contrast to anti-symmetric states to fermions. 

Let us assume three photons labelled as 1,2 and 3 in the same state ip. The three 
photon state can then be written as 

|*> = |Vi,V2,^>- (2) 

In the Fock state representation, this state is denoted by |3). Subtracting a photon from 
this three photon state is due to subtracting either photon 1, 2 or 3. This means that 
the subtraction operation is represented by a symmetric operator / = f\ + f 2 + fo where 
subtracting photon i is represented by fi. Applying / onto |\1/), we obtain 

/>) = |V>2,V>3> + |Vl,'03> + |V'l, V>2> (3) 
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which is not normalised. After normalisation 



/>> = Vs 



-=(1^2,^3) + HlM + \*PlM) 



(4) 



The state described in the square bracket is the state of two photons, i.e. |2) in the Fock 
representation. We thus find that by the subtraction of a photon, the state becomes 

/>> = V3|2>, (5) 

where a/3 is the coefficient required by the annihilation operation in ([1]) when n = 3 
for the state in Eq. (J2J). It is clear that the coefficient v3 appears naturally from the 
argument of the symmetric nature of bosonic particles. 

The coefficient of y/n appears in Eq.flTJ, in fact, tells us that 'the probability of 
a transition in which a boson is absorbed ... is proportional to the number of bosons 
originally in the state 1 as Dirac writes [3j. Similarly we can prove that the creation of a 
photon in state \n) gets the coefficient \fn + 1 as in Eq.flTJ purely using the symmetric 
nature of the bosonic particles. 

On the other hand, we know that in harmonic oscillator, whose dynamics has a 
periodic nature like a wave, the lowering and raising operators, respectively to go down 
and up the energy ladder also have the same mathematical structure as in Eq. flTJ (see 
for example, Chapter 6 of Ref. [1]). Observing that the two are described by the same 
mathematical formulae but looked at from two different points of view, Dirac writes 
this equivalence as 'one of the most fundamental results of quantum mechanics, which 
enables a unification of the wave and corpuscular theories of light to be effected. 1 [3]. 

We can easily derive the commutation relation between the annihilation and 
creation operators using Eq.flTJ: a^a ^ aa) but [a, a)] = 1. In this tutorial, we 
show how to realise the operation of annihilation and creation of a photon. We then 
show how to engineer a quantum state using the properties of annihilation and creation 
operations. 

2. Brief review of quantum optics 

We start with a brief review in order to make this tutorial as self-contained as possible. 
The annihilation and creation operators are not Hermitian, i.e. a) ^ a, as is clear in 
Eq.flTJ, which means that these operators are not measurable. Thus we define quadrature 
operators as linear combinations of the two operators 

q = a-\-a) , p = i(a) — a) (6) 

which are Hermitian and measurable. The m-th order moment of p is 

P m = * m ^ Vp l-o ■ (7) 

Similarly, we can find q m . By defining the characteristic function for a field with its 
density operator p as 

C{u,v) = Tr[pe- ivp+iu «} (8) 
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we can calculate the m-th order moments of p and q at once: 

am flra 

(P m ) = i m -^C(u,v) \ u=v=0 , (^) = (-i)^—- C (u,v)\ u=v=0 . (9) 

The characteristic function C(u,v) is called the Weyl characteristic function. 

In quantum optics, the exponential operator in Eq. (JSj) is called as the displacement 
operator, -D(£): 

-D(0 = exp(— ivp + iuq) = exp(£a* — £*a) (10) 

where the complex number £ = (y + iu). In phase space composed of two conjugate 
variables q and p, D(£) displaces a state by (Re[£], Im[£]) hence the name [5]: 

£>t(£)a£(£) = a + £. (11) 

In probability theory, the characteristic function is Fourier-transformed to give a 
probability density function to tell us about the probability of finding a system at the 
certain value of p and q. However, because of the uncertainty principle in quantum 
mechanics, it is prohibited to simultaneously allocate certain measurement values to 
conjugate variables. Thus, in quantum mechanics, the probability density cannot be 
defined in phase space. By the Fourier transformation of the Weyl characteristic 
function (jHJ), we get a function which is like a probability density function so to call 
it a quasiprobability function [6]: 

W(a) = ^r [ C(u,v)e- 2iarU+2ia * v dudv 
it 2 J 

= ^ J C(0e^- a ^d% (12) 

where subscripts r and j refer to the real and imaginary parts (this convention is used 
throughout the paper). The quasiprobability function was first introduced by Wigner [7] 
as an effort to define a probability-like function in phase space. He derived it under the 
condition that the marginal probabilities / W(p, q)dq and J W(p, q)dq of the probability- 
like function W(p, q) has to be the true probability functions for the respective variables 
p and q. The Wigner function has a one-to-one correspondence to the density operator 
and identifies a physical state uniquely. In practical reasons, the Wigner function is 
very often used to represent the state of a physical system. 

The vacuum state, which is the energy ground state, of a field is denoted by |0). 
Using Eq. (TT2l . we can find its Wigner function in the Gaussian function. By applying 
the displacement operator D(/3) on |0), the peak of the Wigner function is moved by (3 
and we obtain 

W m (a) = ^exp[-2|a - (3\ 2 ] = W ]0) (a - f3), (13) 

where H / |o}(«) is the Wigner function for the vacuum. Glauber [8] introduced the 
displaced vacuum to analyse the higher-order coherence and called such the state as 
the coherent state: \(3) = D(/3)\0). An ideal laser is assumed to be in a coherent 
state even though there is some dispute because the phase of a laser field is completely 
unknown (See [H]). 
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As an example, let us consider to measure the quadrature phase of a superposition 
of two coherent states with tc phase difference: 



'*>= \I^TL \ WS^ + e^ cfc |-/3)), (14) 

n [l + cos(j) sch e 2 p ) 

where (3 is assumed to be real. This state has been studied in conjunction with the 
Schrodinger paradox [22] (See [TU] for a review) and will be discussed further in the 
latter Sections. As an example, when sc /j = 0, the Wigner function of the coherent 
superposition state is calculated as 

W (a) = - l ._,. [e- 2 l a ^l 2 + e" 2 ^! 2 + 2e" 2 H 2 cos^a,)} . (15) 
7r(l + e Ap ) L J 

It is clear that the Wigner function can be negative at some points of the phase 
space, which confirms that the Wigner function is not a probability function. In fact, 
the negativity in phase space is considered to be a signature of nonclassicality for a 
state. There have been attempts to observe the negativity in experiment, which will be 
discussed later. 

Another important state of light to mention is a squeezed state. By a x^ nonlinear 
interaction, a photon in the input pump field is converted into two photons conserving 
the energy and momentum. Arranging the phase-matching condition, the two photons 
can be emitted into the same mode a, in which case the interaction Hamiltonian is 
H s \ = h{\\a? + \\a) ). The strength of the pump laser and the nonlinear coefficient 
determine the coupling parameter Ai. On the other hand, if the two daughter photons 
are emitted into two distinctive modes a and b, the interaction Hamiltonian becomes 
H S 2 = hi^X^ab + \2dJtf). The dynamics of the field in modes a and b can be obtained 
by applying the evolution operator depending on the Hamiltonian. Assuming that the 
field is initially prepared with nothing, i.e. in the vacuum state |0), at the interaction 
time t it transforms into 

|«S 1 }= ( 5 1 (C)|0) = e-^^|0) ; ( = 2i\ l t (16) 

in the case of single-mode emission. The evolution operator S\(() is called the single- 
mode squeezing operator and the state in Eq. (JIB"]) is called the squeezed state (or, more 
precisely, squeezed vacuum). There may be other types of squeezed states by applying 
the squeezing operator to various initial states. In this paper, however, we will be 
interested mainly in the squeezed vacuum. After the decomposition of the squeezing 
operator, we can calculate the precise form of the squeezed vacuum in the Fock state 
basis (See the derivation in [TT] and Appendix 5 in [T2]): 



IS,) = v ^chcf^^2(-itanhC) n |2n), (17) 

n=0 

where ( is considered to be real for simplicity. It is immediately recognised that there 
are only even numbers of photons to be realised in this state, which reflects the nature 
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of twin photon generation by the \ nonlinear interaction. The unitary transformation 
of the bosonic operator is 

Si(C)aSl(C) = a cosh C + a+sinhC . (18) 

If two photons are emitted into two distinctive modes, the vacuum state evolves 



into 



\S 2 ) =,S 2 (C)|0,0) =e-^ s2t/?i |0,0) = sechC^(-tanhCrhri), 



(19) 



n=0 



where the evolution operator is now the two-mode squeezing operator S 2 (C) with the 
squeezing parameter £ = i\ 2 t. Again, the decomposition of the two-mode squeezing 
operator S 2 (C) has been used and ( has been assumed real for simplicity. We observe in 
the two-mode squeezed state (Tl9l that if there are n photons in mode a, we can say that 
there are also n photons in mode b without having to measure it. This 'deterministic 
nature' of the quantum state is an important ingredient of quantum correlation, so- 
called entanglement, the detail of which will be further discussed later. The unitary 
transformations of bosonic operators are 

S 2 (()aSl(() = acoshC+^sinhC ; S 2 (()bSl(() = ocosh(+tf sinhC.(20) 

The generation of an entangled state from a single-mode squeezed vacuum has been 
discussed in [13] (See a review article [14]). 

Using Eq. ( 1TT1) in the definition of the Weyl characteristic function (jHJ), it is clear 
that the squeezing operator transforms the coordinates of phase space, expanding one 
axis at the expense of contracting the other. The Wigner function for the single-mode 
squeezed state is calculated to be 

W s \ 0) {a) = W l0) {a r e (: ,a i e- (: ). (21) 

The Wigner function thus appears squeezed in phase space hence the name squeezed 
state. In fact, the unitary transformations (TTBl and (1201) are the linear transformations 
of the set of initial bosonic operators to that of the output bosonic operators: in other 
words, we can write the transformation in the following general form with the actual 
form of the transformation matrix T to be determined accordingly, 

I a \ 1 a \ 



£1,2(0 



a' 
b 

V Sty 



sUO 



a' 

b 

V 6ty 



(22) 



Huang and Agarwal have shown that the initial Gaussian Wigner function remains 
Gaussian after a linear transformation of bosonic operators of fields |15j . They then 
showed that the action of a beam splitter, squeezer, linear amplifier are all described by 
a linear transformation. With the reasons of mathematical handedness and experimental 
relevance, a Gaussian state in single-mode and multimode fields have been extensively 
studied. 
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Throughout the paper, a homodyne detector is used to measure the properties of 
the light field. The quantum theory of homodyne detection was originated by Yuen and 
Shapiro [T7j. In order to measure quadrature variables, the homodyne detector makes 
use of interference between a reference field, called a local oscillator, and the field to 
probe, by superposing them at a beam splitter. Let us consider a lossless beam splitter 
of two input ports of modes a and b and two output ports of a out and b out . Using 
the SU(2) symmetry, Campos, Saleh and Teich [16] have developed the beam splitter 
operator 

B{e,(j>) = exp [|(a t 6e^ - oSV**) j (23) 

with the reflectivity r = sin | and the transmittivity t = cos |. It is important to note 
that the beam splitter operator is unitary to conserve energy between input and output 
fields. Due to the unitary nature, in fact, two identical photons respectively put into 
two input ports interfere at the beam splitter to exit together at an output port. This 
is the basic idea of the so-called Hong-Ou-Mandel interferometer which appears often 
in quantum technology to prove the identity of two photons [18]. In this paper, we fix 
the phase difference given by the beam splitter to <fi = 0, for convenience. The unitary 
transformation by a beam splitter is then 

u (24) 

which is again a linear transformation. 

Photon numbers are measured at the beam splitter output ports and the difference 
between them is 

h b -n a = (t 2 - r 2 )(a ] a - b ] b) + 2tr(a ] b + ab ] ), (25) 

where the photon number operators n a = a' out a out and n& = b out b out . For a 50:50 beam 
splitter with r = t = l/v2, the photon number difference is 

n a — n b — &b + ah* . (26) 

The homodyne detector with the 50:50 beam splitter is specially called as the balanced 
homodyne detector. Now, we take a strong laser field as the local oscillator which means 
the local oscillator can be considered as a classical field to replace the bosonic operators 
with its dimensionless amplitude El and phase <pL'- b ~ |£|e^ L and b* ~ \El\g~ 1< ^ l ■ This 
approximation brings 

n b -n a = \£ L \ (aV^ + ae~^) = JV^. (27) 

It is easily seen that the photon number difference is the scaled quadrature operator. 
When (pi = 0, iVo = \Sl\(J- When <pL = 7r/2, Nn/2 — \£l\P- Now, it is clear that 
quadrature operators can be measured by the homodyne detector. In fact, by changing 
4>l, we can measure not only two conjugate operators p and q but also the phase 
value along any rotated axis in phase space. By defining the quadrature operator 
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q^ = (ae'^ + ale^), we have the correspondence of the measured value to the quadrature 
variable: N^ L = \8l\q^>- 

According to the probability theory, the probability of photon number difference is 
the Fourier transformation: 

1 



P(N^) = - J dkC H (k)e-"»** (28) 



where 



Ow = Tr 



~ik\£ L \(ae- i 'l'L+ate i 't'L) 
fJK 



(29) 



Comparing Eq. (l2"8|) with Eg. (1121) . we can find that P{N^ L ) is proportional to the 
marginal Wigner function. The probability of photon number difference is the marginal 
Wigner function whose axis is determined by the phase of the local oscillator. With 
this observation, Risken and Vogel proposed to reconstruct the Wigner function using 
the homodyne measurement outcomes by rotating the local oscillator phase [19]. This 
is the so-called optical tomography. 

The homodyne detector is also known for its high efficiency. Without approximating 
the local oscillator as a classical field, Braunstein assessed the reliability of the homodyne 
measurement as a measurement of quadrature phase, depending on the intensity of the 
local oscillator [21]. He found that the strength of the local oscillator is not sufficient 
to ensure that a homodyne detector acts like an ideal detector of quadrature phase. 

Before we finish this section, let us go back to the Wigner function. As it is 
not possible to define a probability density function in phase space, we had to define 
a quasiprobability function which is a probability like function. As it is not the 
true probability function, it is not surprising to have more than one quasiprobability 
function [21] . Using the Campbell-Baker-Hausdorf theorerdjj the displacement operator 
in Eq. (fit)]) can be written as 

D(0 = e ^e-«* a e-l«l 2 / 2 = e ^ & e^e^ 2 (30) 

The Weyl characteristic function Eq. (jSJ) can then be written as 

C(£) = Tr [e^e-^p] e^ 2 ' 2 = Tr [e~^e^ p] e^' 2 (31) 

from which we can see that the expectation value of the normally ordered operators 

<a r <n=(!)"(-£)"ctt).»w (32) 

and that of the antinormally ordered operators is 

<« r Kir(-^)" cKKiitf <33) 

This shows that by Fourier transformation of C(^)e ± '^' / 2 , we will have probability- 
like functions related to the expectation values of normally and antinormally 

\ If [A, [A,B]} = [B,[A,B]} = 0, then exp(0ii + 6 2 B) = exp(9iA)exp(9 2 B)exp(-[6iA,e 2 B}/2).-See, 
for example, p. 42-44 of [12] for a proof. 
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ordered operators. In general, we define s-parameterised characteristic function and 
quasiprobability functions as 

C^s) = C(0e^ 2 = Tim,s)p} ; D(£, s) = D{^)e^ (34) 

where — 1 < s < 1 and 

1 f°° 
P(a, s) = — d 2 £C(£,p) exp«* - «*£)■ (35) 

" J — oo 

It is straightforward to show that P(a, s) satisfies one of the conditions for it to be a 
probability function: J_ P(oe, s)d 2 a = 1 for any s. However, P(a, s) does not have to 
be positive according to the definition. 

When s — 0, P(a, s) becomes the Wigner function. When s — 1, P(a, s) refers to 
the quasiprobability for normally ordered operators and called as the Glauber-Sudarshan 
P function to honour its inventors [Hi [25] (Throughout the paper, the Glauber Sudarshan 
P function is referred by the P function for simplicity). It has been shown that the P 
function is related to the density operator through the following equation 

poo 

p= P(a)\a)(a\d 2 a. (36) 



It is known that any moments of photon number operators should be normally 
ordered [26], for example 

n = (n) ; n 2 = (: h 2 :) = (a i a ! aa) (37) 

where : : stands for normal-ordering. The photon number variance is then 

(: An 2 :) = (: h 2 :) - (h) 2 (38) 

If the P function is positive then P(oe) can be considered as a probability function and 
the normally ordered variance will be positive. However, if P function is not a well- 
behaved positive function, the normally ordered variance does not have to be positive. 
The non-existence of the positive well-behaved P function is thus considered to be a 
sign of nonclassicality for a given field. 

Now, we have seen two criteria of nonclassicality: one is based on negativity of 
the Wigner function and the other is based on the non-existence of the well-behaved 
P function. These two conditions do not necessarily coincide so there seems to be a 
subtle problem to be settled. In fact, as we can show later, the positivity of P function 
guarantees the positivity of the Wigner function but the converse is not necessarily true. 

We can obtain P(a, s') from another quasiprobability function P(a, s) for s' < 

P(a, s>) = — *—, [^ d 2 (3e W f ^-P} 2 ) P(/5, s ), (39) 

v(s - s') J_ 00 V s ~ s J 

which appears as a convolution of P(a, s) with a Gaussian function and its effect is 
to bring the peaks of P(a, s') broader than those of P(a, s) (In fact, this broadening 
reminds us of a decoherence process of the field and the detail has been studied using 
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the Gaussian convolution theory [28]). The s-parameterised quasiprobability function 
( 135]) together with (|34p can be written as 

P(a, s) = Tr[fr (a) pD (a) f(s)] (40) 

where T(s) — -j J d 2 £.D(£, s). It can be shown for s < 1 that [12] 

f(s) = l : exp ( --?—tf a ) : (41) 

77(1 — Sj \ 1 — S / 

Then expanding the exponential function and using the completeness J2T=o \k)(k\ = 1> 
we find that 

: exp I — a' a 

= 1 - a T " -i — I — — i » T 7r -77 ( — — i «' <'" i 

oo fc 

= ££(-7-) i;)i*><*l = E(^)l*X*l- («) 

fc=0 n=0 V / \ J k=0 

Substituting this into Eq. (J4"D1) . the s-parameterised quasiprobability function can be 
written as 

P{^s) = ——-Y J [-^ l ) (k\D\a)pb{a)\k). (43) 

which allows us to have a connection between the quasiprobability functions and the 

density matrix elements. In fact, through this connection, we can also reconstruct 

a quasiprobability function after the density matrix is reconstructed from homodyne 

measurements (See [22] for details). 

When s = 0, P(oe, s) in Eq. (143 j) becomes the Wigner function: 

o °° 
W(a) = - Y(-l) k {k\D^(a)pD(a)\k) (44) 

77 *■ — ' 
fc=0 

The Wigner function at a of phase space is the probability of having an even number 
of photons less that of having an odd number of photons after displacing the field by 
a. The dependence of the Wigner function on the even and odd parities is important 
for its reconstruction and applications. When s = —1, P(a,s) is called the Q function, 
which is the quasiprobability function related to the expectation values of antinormally 
ordered operators. From Eq. (1431) . we can find that Q(a) is the overlap between the 
field and the coherent state \a): 

Q{a) = -(a\p\a). (45) 

77 

Thus the Q function is positive at any point of phase space. At one end s = — 1, the 
s-parameterised quasiprobability function is positive well-behaved for any state. As 
s grows, the peaks of the function become sharper. At the other end s = 1, the s- 
parameterised quasiprobability function may not even exist depending on the state of 
the field. 
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For quantum state reconstruction, the homodyne measurement with dark count 
and other sources of inefficiency have to be assessed. The optimum number of 
measurement angles also has to be assessed. After the data is collected, the quantum 
state reconstruction has to be considered with the algorithm of maximum likelihood [22]. 

3. Photon subtraction and addition using a cavity or a beam splitter 

While photon addition and subtraction have been discussed for long, neither of the 
operations were experimentally demonstrated till recently. Wenger, Tualle-Brouri and 
Grangier [29] were able to subtract a single photon from a travelling field using a 
beam splitter and a photodetector. In the same year, a single-photon addition was 
experimentally performed by Zavatta, Viciani and Bellini [30]. The combination of the 
subtraction and addition has been successfully demonstrated in [31]. 

In order to illustrate how to add and subtract a photon, let us consider a perfect 
single-mode cavity where a state of a field has been prepared (see Fig. [I] (a)). We 
then send a two-level atom of its excited |e) and ground \g) states through the cavity. 
Assuming resonant interaction between the atom and the cavity field, the Hamiltonian 
in the interaction picture is written as 

Hj c = h~\j C (a)d-- + aa + ), (46) 

where Xj C is the coupling constant, <j+ = \e)(g\ and <7_ = \g)(e\. This fully-solvable 
Hamiltonian has been widely studied under the name of Jaynes-Cummings model (See 
[32] for a review). Let us assume that the atom is injected into the cavity in its excited 
state then it is measured in its ground state after interaction with the cavity field at the 
exit of the cavity. This is described by 

A jc = (g\e-^\e) = _ ^(W^ a t. (47) 

where the orthonormality (g\e) = has been used. When 3< j- n j ^ ^ ^ e a t m 
will add a photon to the cavity field without any extra operation. The condition is 
satisfied when Xj C is small and the photon number inside the cavity is small. Similarly, 
if the atom is prepared in its ground state and measured in its excited state after the 
cavity field interaction, a photon is subtracted from the cavity under the same condition. 
The operation is described by Sj c = —ia 3 % - ■ In fact, a single photon state has 
been generated experimentally by sending an atom to an empty cavity (prepared in the 
vacuum |0)) [33]. Note that any constant value in front of a or a) is neglected because 
that has to disappear when the final normalisation is taken place. Throughout the 
paper, we will ignore those constant values. 

However, according to the input-output theory, the photon subtraction and addition 
assisted by a cavity does not seem to be an optimum solution even though a cavity may 
be a candidate for other types of light field engineering. Let us thus take a travelling field 
in mode a to a beam splitter as shown in Fig. [1] (b), where nothing has been put into the 
other input port of mode h ('nothing' means the vacuum state). If the injected vacuum 
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(a) 




Cavity 



(b) 




inefficient PD 



Figure 1. (a): Operations on a field by passing a two-level atom through a cavity. The 
atom is initially prepared in its ground state and then measured in its excited state 
at the exit in order to subtract a photon (filled dots in the level scheme). The atom 
is initially prepared in its excited state and then measured in its ground state at the 
exit in order to add a photon (empty dots in the level scheme); (b): Operations on a 
field by using a high-transmittivity beam splitter. A photon is subtracted conditioned 
on the vacuum input to take away a single photon and detected as a single photon 
state at the output. If a single-photon state is input to the beam splitter and leaves 
it as a vacuum, a single photon has been added to the field input to the other input 
port of the beam splitter. In the inset, the equivalence of an imperfect photodetection 
with its efficiency g is depicted using a perfect photodetector with a beam splitter of 
transmittivity g in front. BS: beam splitter. PD: photodetector 



takes away a single photon, the action of a photon subtraction has been performed. In 
order to analyse this action, let us consider the generators J + = ab> and J_ = a)b of the 
beam splitter operator B(9, <fi = 0) in Eq. ( 1231) . Together with J 3 = \(jftb — a^a), J + and 
J_ form a group thus we can decomposq§| the beam splitter operator as (For example, 

§ The Schwinger's relation |36] is to cast the two-dimensional harmonic oscillator in terms of an 
angular momentum system normalised to h. The relation is based on the fact that J\ — oW+ + ■'— )) 
Ji = j[(J+ ~ J-) an d ^3 have the same commutation relations as the angular momentum operators. 
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see pages 45-46 in [12]) 



B($,0) = e " tan 2 J + e -ln(cos 2 5 )J 3 e tan 5 J__ /^g\ 

As input mode b is in the vacuum, 

B(^ ) 0)|0) fe = e- tan ^ J +t ata |0) b (49) 

where t is the transmittivity of the beam splitter as in Eq. ( |23l) . This is further reduced 
if we project it onto a single-photon state. Thus the action of single-photon annihilation 
is summarised as 

S bs = b (l \B(9, 0)\0) b = - tan -a£ ata ~ a£ ata . (50) 

It is clear that the projection gives an impact of subtracting a photon from the field in 
mode a but in order to leave only the annihilation operator, the following condition has 
to be satisfied 

t ata « 1, (51) 

which can be compared with the condition leading to the annihilation operation 
using the Jaynes-Cummings model. The factor — tan | in Eq. ( 1501) disappears during 
normalisation. According to the definition of the beam splitter operator, the parameter 
9 determines reflectivity and transmittivity, and 9 small means the reflectivity small. 

If, instead of the vacuum, a single-photon state is injected to mode b and nothing 
is measured at the output, we know that one photon has been added and the overall 
operation is 

A bs = b (0\B(9,0)\l) b = tan V^ ~ £ a V. (52) 

Under the condition ( 15TT) . the projection A bs becomes to simulate a) . However, this 
scheme is demanding because, the single photon generation on demand is challenging 
and the detection of the vacuum accompanies higher noise than the detection of photons. 
In the photon subtraction scheme using the beam splitter, the most important 
experimental challenge is the single-photon detection. There are many issues in the 
inefficient detection including 1) the low detection efficiency of a photodetector and 2) 
the detection of a photon in a wrong frequency band. An inefficient photodetector with 
its detection efficiency n can be simulated by a combination of a perfect photodetector 
with a beam splitter in front (See Fig. CD (b)), where the vacuum (mode v) is assumed to 
be injected to the unused input and the transmittivity of the beam splitter is the same 
as the detection efficiency n. Assuming the input field to subtract a photon is p a , the 
output field in mode a is 

Ps = Tr v ( b (l\B bv B abPa ® (10,0)^(0,01)5^11),) (53) 

Using the expansion of the beam splitter operator (!2"3l under the condition (l5~Tj) . 



P'™\2J ^aP^+U) r, 2 (l-n 2 )a 2 p a a) 2 . (54) 

As the operators form a complete group, it is clear that the decomposition is possible even though the 
actual decomposition is not necessarily straightforward [37) . 
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The first term represents the annihilation of a single photon while the second term is to 
annihilate two photons. For the second term to be negligible, the reflectivity (9) of the 
beam splitter has to be small and the detection efficiency has to be large. 

Currently, the efficiency of single-photon level detection is very low, even though it 
depends on the frequency of the field to be probed. On the other hand, an avalanche 
photodiode can be used to detect photons as it is highly sensitive to a photon. It is, 
however, saturated by a single photon, which means that a photodiode is like a on-off 
detector to be 'on' if there are photons and 'off' if there is no photon. It cannot tell 
how many photons there are. Mathematically, 'on' is to project the field onto 

oo 

G = X>)<n| = l-|0><0|. (55) 

Instead of a single-photon detector, the on-off detector is used in recent photon 
subtraction experiments [29], [Ml G22, [30] . In order to see how it works in the annihilation 
of a photon, let us replace the single-photon detector with G in Eq. f[50|) . For the density 
operator for the input field in mode a is p a , the output conditioned on an 'on' event is 

Tr b [G b B(9, 0)p a ® (|0) 6 (0|) B\6, 0)] = Tr b [B(9, 0)p a ® (|0) 6 (0|) &{0, 0))] 

b (0\B(9,0)p a ®(\0) b (0\)B\9,0)m 

where 

Tr b [B{9, 0)p a ® (|0) 6 (0|) B\0, 0))] « £ ata p a £ ata + tan 2 °-a £ ata /3 a t ata a 1 " (57) 

under the condition (j5lj) and 

b (0\B(9, 0) Pa ® (|0) 6 (0|) B\9, 0)|0} 6 = t ata p a t ata (58) 

Substituting the above calculations into Eq. ( |56l) . we find that 

Tr b [G b B(9, 0)p a ® (|0) 6 (0|) B\9, 0)] « tan 2 °-a t ata p a t ata a f . (59) 

Like in Eq. (1541) for the consideration of photodetector inefficiency, the use of an on-off 
detector inevitably results in the conditioned state to be mixed but under the usual 
condition, the resultant state is nearly a single-photon annihilated state. 

We note that all the schemes to add and subtract a photon are conditional based 
on projective measurements. A photon is added and subtracted conditioned on the 
initial preparation and final measurement outcome. In order to simulate annihilation or 
creation, we have found that to put or to take away extra unit of energy is not enough 
but the interaction with the process has to be weak. Otherwise, the field state gets extra 
kick due to involvement of higher-order terms in the Taylor expansion of the interaction 
operator. 

3.1. Quantum scissors 

Before we finish this section, we would like to briefly review a quantum engineering 
scheme based on the projective measurement like photon annihilation and creation 
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introduced above. Earlier, Pegg, Phillips and Barnett [3H] pioneered a scheme called 
the quantum scissors for state truncation. The quantum scissors are composed of two 
concatenated beam splitters and photodetectors. The output mode b of the first beam 
splitter is arranged to be an input to the second beam splitter. A single photon is input 
to mode a of the first beam splitter while no photon is input to mode b and a field of 
state \ipsc) to input mode c, where c is to denote an input mode of the second beam 
splitter. If the state is \ip sc ) = 7o|0) +7i|l) +7 2 |2) + ■ ■ ■ then the final state at the three 
outputs is 

^|l)„|0) 6 |0) c - ||0) a (|l) 6 |0) c - |0) b |l) c ) + ||l) a (|0) 6 |l) c - |l) 6 |0) c ) 

-y|0)a(|0) b |2) c -|2) fe |0) c ) + ---.(60) 

Now, conditioned on the measurement outcome and 1 at modes b and c respectively, 
we effectively select the state in mode a, as 

^(yl°)a + y|i>.) ( 61 ) 

where M is the normalisation factor. We see that the input |l) a has turned into a 
superposition state whose composition depends on the input state in mode c. In fact, 
taking \ipsc) as a coherent state we have a fair amount of freedom to choose 70 and 71. 

4. Quantum-state engineering by photon addition and subtraction 

Dakna and coworkers [JT] proposed a scheme to generate a travelling field of an arbitrary 
superposition of Fock states. Their scheme is based on the sequence of displacement 
operation D(a) and addition of a single photon. As a simple example, let us consider to 
generate a superposition of |0) and |1), which may be compared with a quantum scissors 
scheme to generate such the superposition as in Section I3TT1 The arbitrary superposition 
state co|0) + ci|l) with arbitrary complex numbers Co and C\, can be written as 

ci(^|0> + |l>). (62) 

Using the creation and displacement operators, state (1621) can be written as 

Cl (V + |) |0> = Cl L>t f^\ a^D h) |0> (63) 

where Eq. ( fTT|) has been used. Thus, sending a vacuum to a sequence of a displacement 
operation, single-photon addition and another displacement operation, we obtain the 
state (j52j). 

Using the equivalence of any superposition state 

N N 

i*> = E c '>> = E7V at ) n i°> (64) 

n=0 n=0 V n - 

to 

|tf > = nJUCot + <)|0) = Il» =1 fr(a n )tfD(a n )\0) (65) 
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Dakna and coworkers [JT] showed that an arbitrary superposition of |0), |1), • • • , \N) can 
be generated by a concatenation of N units of operations composed of displacement, 
single-photon addition and another displacement as shown in Fig. [2j 

The displacement of a quantum state in phase space can be realised [12] using a 
beam splitter of high transmittivity t « 1 and a strong coherent field of its amplitude 
(3. Let us consider a beam splitter operation where one input mode is a coherent field 
of \j3) b : 



B( 



. = 0)|/?) b = e-^ b - ab '>\f3) b ~ e ^ a,p - ap } \(3) h (66) 

where the approximation was done for the condition of a large amplitude /3. The result 
of the operation is the displacement of ~/3 for mode a. More rigorously, we can discuss 
the displacement operation as follows. Any field state of density operator p can be 
written as a weighted sum of coherent states as in Eq. (136]) . Thus by beam splitting a 
field of an arbitrary state in mode a and a coherent field \(3) in mode b, we have 

p„b = B(e,0)p a ®\P) b (P\&{9,0) 

= £(0,0) f d 2 aP(a)\a) a (a\ ® \/3) b (J3\&{e,0) 



d 2 aP(a)\ta + r(3) a (ta + r(3\ ® |t/3 - ra) b (tf3 - ra\. (67) 

The state of mode a regardless of mode b can be obtained by tracing over mode b: 

Tr 6 p ab = / d 2 aP(a)\ta + r(3) a (ta + r(3\ (68) 

which is approximated to the displaced state 

D{r(3)p a D\r(3)= f d 2 aP{a)\a + r(3) a (a + r/3\ (69) 



I0> 



D(a.) 




Displacement 



Figure 2. Diagrammatic sketch of the scheme to generate an arbitrary superposition 



of Fock states Yl 



N 



cl* denotes the operation of a single-photon addition and 



D(ai) (i — 1, ..., N) is the displacement operation, which can be realised by the use of 
a beam splitter of high transmittivity t and a strong coherent field of amplitude (3 as 
shown in the magnifier. The displacement a = f3\/l — T with T = \t\ 2 
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when t«l. As r = sin | rs | for ?a 0, Eqs. (I66p and (1691) are in good agreement. 
However, as can be seen in Eq. (1681) . the displacement operation is never exact as when 
t — 1, r = and a state is never displaced and when £ < 1, state ( 1681) is never the same 
as ([69]). 

While adding a photon is difficult to realise, subtraction is relatively easy (an 
experiment of single-photon addition will be explained later in this paper). Fiurasek, 
Garcia-Patron and Cerf [13] proposed to engineer a quantum state by a sequence of 
subtraction and displacement of a squeezed vacuum. First, a squeezed vacuum is 
displaced before a single photon is subtracted from it. After another displacement 
operation, the state is antisqueezed. All these processes are summarised as follows: 

\*)aut = 5 1 t (C 2 )^(a 2 )a J D(a 1 )5(Ci)|0). (70) 

If a,\ = a 2 and Ci = C2 then 

|^) ou< = 5 1 t (Ci)(a + ai)5(Ci)|0) = -sinhCi|l) + ai|0), (71) 

where the unitary transformation (TT8]) has been used. By an appropriate choice of (1 and 
«i, we can have an arbitrary superposition of |0) and |1) states. It can be proven [13] 
that we can achieve an arbitrary superposition of Fock states, by concatenation of the 
units composed of subtraction and displacement of the squeezed vacuum together with 
the final antisqueezing. The proof follows a similar analysis shown earlier in this Section. 
In the next subsection, we will show how to produce a specific quantum state using a 
squeezed vacuum. 

4.1. Production of a cat state 

The sequence of photon subtraction or addition is a demanding task in experiments. 
As a simple example of quantum-state engineering, let us consider recent experimental 
triumphs to produce linear superpositions of coherent states, by subtracting a photon 
from a squeezed state. It has been a long awaited dream to produce a coherent 
superposition state (TH1) . which was briefly discussed in Section El A coherent state 
is considered to be at the boundary of classical and nonclassical worlds and its 
superposition is accepted to show the quantum nature of Schrodinger's cat parado^ljj]. 
There has been a difficulty in producing such the superposition state in a travelling field 

|| In his argument Schrodinger tried to show how the microscopic world which may be under the laws of 
quantum mechanics is directly coupled to the macroscopic world which is under the laws of conventional 
physics at the time. In this paradox, a cat is placed in a steel box and a poisonous gas is set to be 
released when an atom, connected to the container of the gas, decays. Reading |39j . it is not clear 
that the state described in the paradox is an entangled state between the atomic state and the destiny 
of the cat or it is a simple superposition of the cat's destiny. Had it been the former, the state for 
the cat and atom will be \no — decay) \alive) + \decay)\dead) assuming the equal amplitude of the two 
events of atomic decay and non-decay. Here, normalisation has been neglected. If the paradox is simply 
described by the superposition of a cat state, then \dead) + \alive) will describe the paradox. 
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because of the unavailability of extremely high nonlinearity required 
state is a Poissonian weighted superposition of Fock states [8j 



The coherent 
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thus assuming the amplitude j3 to be real, the coherent superposition state (H 
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(74) 



V^sinh/? 2 ^ v^n + l) 

for 0^ = 7r. It is clear that |\l/ e ) (|^o)) has the non-zero probabilities only for even 
(odd) numbers of photons hence the name the even (odd) coherent state. The Poissonian 
weight is shown in Fig. [3J 

Recalling that the squeezed vacuum has non-zero probabilities of having only the 
even numbers of photons, as shown in Eq. (1171) . we can immediately recognise some 
analogy between the photon number distributions of the coherent superposition state 
and the squeezed vacuumjj. While the weight function in Eq. (TP7J) is a decreasing 
function with the peak at n = 0, that in Eq. (1731 is Poissonian with the peak at (3 2 . If 
we can push the peak of the weight function (1T71) . it seems to be possible to get the two 
distributions closer to each other. By annihilating a photon from the squeezed vacuum, 
the state becomes (without normalisation) 

aS\0) = SS j aS\0) = S{acosh( + aUinh()\0) -> Stf\0) (75) 

% We analyse the similarity of two states based only on their photon number distribution in this part of 
the paper. This analysis is valid when the two states are pure and the weights of Fock-state components 
are real and positive. Although it is dangerous to identify a quantum state under such the restrictive 
constraints, it provides an intuitive picture so to use it in this part of the paper. 
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Figure 3. Photon number distributions for the even coherent states of (a): (3 = 2 
(b):/3 = 3. 
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where the unitarity of the squeezing operator has been used together with the 
transformation properties (IT51) . The argument £ of the squeezing operator has been 
omitted to make the notation simple and the sign — > means that its left hand side 
becomes the same as the right hand side after normalisation. We note that annihilating 
a photon from the squeezed vacuum gives the same result as squeezing a single photon 
state. Furthermore, we will see that annihilating a photon in the squeezed vacuum gives 
the same impact as adding a photon to the squeezed vacuum. By applying a creation 
operator to the squeezed vacuum, we get 

tf§\0) = S'^a+S'lO) = 5(a t coshC + asinhC)|0) -► Sa^O). (76) 

It is surprising that adding a single photon to a squeezed vacuum gives the same 
impact as annihilating a single photon from it. A starting point to understand this 
is the fact that there are only even numbers of photons in the squeezed vacuum. By 
subtracting a photon, we know that there were at least two photons in the field, as 
otherwise it would not have been possible to subtract any photon. By subtracting a 
photon, we now know that there is at least one photon in the field. Thus, subtracting 
a photon, we annihilate a chance of having no photons and results in the same effect as 
adding a photon to the field. Here, we remind ourselves that the process is conditional. 

By subtracting a photon, we have been able to move the peak from zero photon to 
nonzero photons. This is well seen in Fig. HI The position of the peak depends on how 
many photons are annihilated and how much the state is squeezed initially. The reason 
why the peak moves to a higher photon number is because of the coefficient \/n of the 
annihilation operation as shown in Eq. (CQ). This coefficient gives a higher weight to the 
initial higher- number photon component. 

For the photon subtracted squeezed vacuum to simulate the coherent superposition 
state, the weight distribution of its component states should be close to that of 
the coherent superposition state. The higher the squeezing, the flatter the weight 
distribution for the initial squeezed vacuum. Thus by subtracting a photon we have a 
longer and thicker tail in the weight distribution than in the Poissonian distribution for 
the coherent superposition state. By an appropriate choice of the squeezing parameter, 
the single-photon annihilated squeezed state is observed to simulate the coherent 
superposition state [H]. A sequence of single-photon annihilations will gradually 
increase the weights of higher numbers of photons (Noting that a coherent state is 
an eigenstate of an annihilation operator, by annihilating a photon from an even 
coherent state, it turns into an odd coherent superposition state and vice versa [49]. 
However, the photon-subtracted squeezed vacuum changes its nature slightly as it is 
not a perfect coherent superposition state. As the number of subtraction operations 
grows the resultant state digresses from the coherent superposition state.). Now, we 
face a problem to measure how close two states are to each other. One widely accepted 
measure is the fidelity which is simply an overlap |(0|0)| 2 between two states \ip) and 
|0). This measure can be used if at least one of the fields is in a pure state. If the other 
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state is mixed so to be represented only by its density operator p^ then the fidelity is 

F = WM) (77) 

which can be calculated using the overlap between their Wigner functions and Weyl 
characteristic functions: 

T = n I d 2 aW^(a)W (l> (a) 



IJd^c^oc^-O 



(7? 



When the two states are the same, the fidelity is 1 with the perfect overlap while when 
the two states are orthogonal to each other, it is 0. The fidelity, however, is not an 
ultimate measure as it can smear out some useful information. Moreover, in general, we 
do not know how to measure the fidelity experimentally, except a few cases [45]. For a 
physical system which is described in a small Hilbert space, the distance measure can 
be more useful and it has its implication in experimental realisation [46]. 

It was calculated theoretically [U] that the fidelity between the coherent 
superposition state and the single-photon annihilated squeezed vacuum is found to be as 
high as 0.99 when the squeezing parameter ( = 0.43 and the coherent state amplitude 
(3 = 1.16. The production of such the superposition state was successfully observed 
experimentally [35], under the name of a Schrodinger kitten state (While Ref. [35] 
reported the production of the coherent superposition state in a pulsed field, Ref. 
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Figure 4. Photon number distributions for (a): the squeezed vacuum (b): the single- 
photon annihilated squeezed vacuum, and (c): the two-photon annihilated squeezed 
vacuum. The squeezing parameters are £ = 1 for (a) and (c) and £ = 0.5 for (b). 
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demonstrated the production of such the state in a continuous wave [48].) because 
the amplitude of the component state is very small, which is a problem to use it for 
many of the tasks in quantum information processing. In their work, Ourjoumtsev and 
coworkers [50] have used a somewhat different scheme to realise a coherent superposition 
state with a larger amplitude as explained below. 

A squeezed superposition of Fock states may be useful for various applications. 
Going back to the Fiurasek and coworker's scheme [13], a squeezed superposition state is 
generated by omitting the final squeezing. In this respect, the generation of a squeezed 
Schrodinger cat state has recently been demonstrated experimentally [50]. Here, the 
superposition of two squeezed coherent states with amplitude a = ±\^L6 and squeezing 
of 3.5dB, was produced. The scheme starts from the generation of a Fock state \n). The 
Fock state is sent to a 50:50 beam splitter. Then a homodyne detector is placed at one 
output port in order to select a state at the other output, conditioned on the homodyne 
measurement outcome q = 0. 

The state conditionally generated is 

\^sc)a = b(q(0)\B (9 = -|, <f> = 0) |n)„|0) 6 (79) 

where the eigenstate \q(x)) of the quadrature operator q can be written as the 
superposition of Fock states 

1 °° 1 
l^)> = -Tu E -^r=f e - x,2 H n {x)\n). (80) 

A Hermite polynomial H n (x) for argument has a non-zero value only for an even index, 
i.e., H 2n (0) = (-2) n (2n - 1)!! and H 2n+ i = 0. Thus 



oc 



W0» = ^E<-ir^|2»> (8!) 

Using the decomposition of the beam splitter operator (jlHj) . the beam splitter operation 
in Eq. fl79j is 




M >:i : 1 1» -*).!*>*■ m 



B(6 = 7i/2,<p = 0)\n) a \0) b 

n n 

■^) ,=0 

Conditioned on q = for mode b, we see that mode a has only even or odd numbers 
of photons depending on the initial preparation of the Fock state. For example, if n 
is even, the state generated is of even numbers of photons. The weight function of 
the Fock-state components is determined by the binomial distribution and the Hermite 
polynomial. As the peak is shifted by the original excitation of the Fock state, we 
can see a possibility of generating a cat state with a larger amplitude. The squeezing 
of the cat state appears due to the final weight function as the overlap of the Hermite 
polynomial and the binomial distribution. We note that under this scheme, the coherent 
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superposition state, which lies in the infinite-dimensional Hilbert space, is realised in a 
n-dimensional Hilbert space restricted by the initial preparation of the Fock state. 

In the experimental realisation, the bottleneck of this scheme is the initial 
preparation of the Fock state. In 2006, the experimental group in Paris [51] successfully 
produced a two-photon Fock state |2) using an optical parametric amplifier. In this 
experiment, the non-degenerate downconversion converts one pump photon to two 
photons into two different modes as explained in the two-mode squeezing operator in 
Eq. ( Tl9l) . At one output mode of the optical parametric amplifier, a 50:50 beam splitter 
is placed as shown in Fig. [5j Detecting a photon at each of the output mode the input 
field is projected onto 

1 :(6i(2| fe2 (0|+ fel (0| fe2 (2|) (83) 



6i(i|m<i|£ 
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As input mode 62 is served by the vacuum, the coincidence count is to project the 
61 mode onto the two-photon Fock state. Due to the twin-photon nature of the 
downconversion, when the field mode b is projected onto the two-photon Fock state, 
the field mode a is also in the two-photon Fock state. In this way, the experimental 
group in Paris produced |2) state and proved it by the full reconstruction of its Wigner 
function. 

4-2. Photon addition using parametric down converter 

We have considered schemes to add a photon using a cavity or a beam splitter in 
Section [31 In this subsection, we show another way to add a photon to a travelling 
field which has recently been realised experimentally. Agarwal and Tara [52] found 
that adding a definite number of photons to a coherent state, the state turns into a 



Pump 




Click 



Figure 5. An optical downconverter produces twin photons into modes a and b±. 
Conditioned on the projection of mode b\ onto |2), we generate the two-photon state 
1 2) in mode a. The condition is based on counting single photons (clicks) in modes b\ 
and hi at photodetectors (PDs). 
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highly nonclassical state, showing a sub-Poissonian character and quadrature squeezing. 
Lee [53J later realised that if a single photon is added to any field state (pure or mixed), 
the state turns into a nonclassical state for which a positive well-behaved P function 
cannot be imposed. He defines the depth of nonclassicality using the properties of the 
quasiprobability functions. As shown in Eq. (1391) . the quasiprobability functions are 
connected to each other. At one end of the spectrum (s = 1), the P function may not 
even exist as a rational function for a nonclassical field while at the other end of the 
spectrum (s = —1), the Q function is always positive and possesses all the properties 
as a proper probability distribution. When s = 1, Eq. (1391) becomes 

1 f°° f — \a - B\ 2 \ 1 - s' 

P^r) = —J d^expf ' - P| JP((3) ; t=-£-. (84) 

Lee defined the depth of nonclassicality based on how well the quasiprobability function 
behaves like a proper probability distribution. The minimum r which brings P(a, r) to a 
proper probability function is the depth of nonclassicality. This means that through the 
spectrum of s, if a state has its quasiprobability function P(a, s) positive well-behaved 
only for its Q function, the state is said to be most nonclassical. This definition is related 
to the resilience of the nonclassicality to the vacuum noise as its influence appears in 
the form of Gaussian convolution (See [28]). Lee proves that for any state of density 
operator p = X]^°m=o/ 9 ( n ' m )l n )( m l' ^ P(0;0) = then the nonclassical depth of the 
state is 1. According to this theorem, even a very chaotic thermal field at a very high 
temperature can become nonclassical only by adding a single photon to it. Of course, 
how useful it can be for any kind of quantum processing is another matter. 

Even though the addition of a single photon (or a definite number of photons) 
guarantees nonclassical behaviour, it had to wait for long before such the operation 
was first realised experimentally by Zavatta, Viciani and Bellini [30]. Realising that 
adding a photon conditioned on 'no ' detection of a photon using a beam splitter and a 
single photon input is experimentally demanding, they add a photon using a parametric 
downconversion process. If a photon is detected in mode b (see Fig. [6]), we know that 
its twin photon is in mode a. If \ip) is an input to mode a, the output field conditioned 
on a photodetection in mode b is Afar\ip) with normalisation Af '. By adding a photon to 
an input coherent state, Zavatta and coworkers claim the observation of the transition 
from a classical to a quantum state. They show this by the full reconstruction of the 
Wigner functions. They also show the nonclassicality of a photon-added thermal state 
in another experiment [54"] . 

The parametric downconversion is described by the two-mode squeezing operator 
f fT9|) . which can be decomposed into (See Appendix 5 in Ref. [12] ) 

q ( /-\ „— a'b* ta.nh(^ — (a* a+W ) \n(cosh £) featanh^ f'&Ei'l 

Detecting a photon in mode b at the output while its input is vacuum, is equivalent to 
the addition of a photon in mode a: 

6<1|&(C)|0> 6 = -^7at(coshC)- ata (86) 

cosh Q 
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where the constant term should be ignored as it will disappear during the normalisation. 
When the squeezing is small, in other words, when the efficiency of the parametric 
downconversion is small, and the input photon number is not large, the term (cosh£)~ a a 
is negligible. Thus the process becomes to simulate the photon creation operation at. 

4-3. A sequence of photon addition and subtraction 

Photon addition and subtraction are now experimentally plausible schemes and in this 
section, we combine the two schemes. Quantum mechanically, photon addition and 
subtraction correspond to bosonic creation and annihilation operations, respectively. 
As the two operators do not commute, the two different sequences of addition and 
subtraction will bear different states. Here, we show how significant the differences are. 
By annihilating and then creating a photon to a coherent state \/3), the state 
becomes 



W) = y/77JFtfa\/3) ; M™ = l/n c (l + n c ) (87) 

where n c = \/3\ 2 is the mean photon number of the coherent field. The photon number 
distribution for the resultant state is 



P* s (n)=N: s nP c {n) ; P c (n) 



n! ' 



which obviously has null zero-photon probability. In order to quantify the uncertainty 
of the photon number in comparison to the Poissonian level, we define the Q- 
parameter [55] . Q = (An) 2 — n. For Q < 1, the field is said to manifest sub- Poissonian 
statistics while for Q > 1, the field is super- Poissonian. Considering the fact that a 
coherent state is an eigenstate of a, i.e. d\j3) = (3\f3) and the photon number operator 
n = a) a = aat — 1, we find the Q-parameter for the field ( 1H71) as Q = — 1 — (n/(n + l)) 2 . 
The field ( 187)) thus is sub-Poissonian. 



Downconverter 




PD 
Click 



Figure 6. Zavatta, Viciani and Bellini [30] successfully added a photon into a coherent 
field using a downconverter. Conditioned on photodetection in mode b a photon is 
added to the input field \i/j}- 
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If we add a photon before subtracting one for an initial coherent state, the state is 



IVT) = VW^IP) ; N s c a = l/{hl + M c +l). (89) 

Its photon number distribution is P^ a (n) = (n + l) 2 P c (n)/(n 2 + 3n + l), which is different 
from the initial P c {n). We find that the Q parameter is negative for any value of /3. 

By subtracting then adding a photon to a thermal field, the density operator 
Pb = Z~2n B /(l + n B ) n+l \n){n\ of the initial field becomes 

p™ = N^a)ap B a)a ; A^ s = l/(n B (2n B + 1)). (90) 

The average photon number and the variance for the resultant state are 

= 2nB(3flB + 2) 2 = l2 _ % + u _ b + 

2n B + 1 
which leads the Q parameter smaller than zero when the average photon number n B of 
the initial thermal state is less than around 0.6. 

The fact that subtracting and adding a photon brings a thermal field to a sub- 
Poissonian state, may inspire us to conjecture that subtracting a photon can narrow 
down the width of the photon number distribution. By subtracting £ number of photons 
from the thermal field, the Q parameter becomes n 2 (£ + 1) which is always positive. In 
fact, not the subtracting but the process of adding a photon turns the thermal field into 
a sub-Poissonian field as the resultant state has the Q parameter equals 2n 2 — 1. 

For the reverse process of subtracting after adding a photon, the state becomes 

p™ = U^aa) p B aa) ; N s B a = l/((n + l)(2n + 1)). (92) 

The average photon number in this case is smaller than that for p a B by one. Even though 
the photon number distribution is not the same as the original one, the resultant state 
is always super-Poissonian with Q > 0. 

In order to see the negativity of the Wigner function, we use a theorem [27] that if 
the Q function becomes zero at any point of phase space, the Wigner function should 
have negative values. Using this, we can easily check if a state is nonclassical because 
the Q function, which is defined in Eq. (HoT) . is much easier to calculate than the Wigner 
function. 

For IV^)) the Q function is found to be 

Kf as 
Q™( a ) = l^|/5| 2 |a|V l/3 - a|2 (93) 

7T 

which is zero at the origin of phase space. Thus, the Wigner function should manifest 
negativity and the resultant state is nonclassical. In fact, we can easily prove that any 
field becomes nonclassical by subtracting and adding a photon. The density operator of 
a classical state can be written as p c i = J P c i(P)\P){j3\d 2 j3, where P c i((y) is a positive well- 
behaved function. The Q function for the field after subtracting and adding a photon is 
the weighted sum of (193]) with a proper normalisation. Thus the Wigner function should 
be negative at some points of phase space as the Q function is zero at the origin of phase 
space. This can be understood easily as the state after photon subtraction and addition 
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is never empty and according to Lee [53] this field of zero null-photon probability is 
nonclassical. 

For \ipc a )i the Q function is 

Q?(a) = -^\a*(3 + 1|V I/3 ~ Q|2 (94) 

TV 

which surely becomes zero when a* (3 + 1 = so the field is nonclassical. However, 
we cannot conclude that any classical state becomes nonclassical by adding and then 
subtracting a photon. It is because 

Q™(a) oc / d 2 (3\a*(3 + l\ 2 e~^- a \ 2 P d {(3) (95) 



does not seem to be zero at any point of phase space. Here, the integrand is a positive 
function and a cannot make the integrand zero regardless of (3. The Wigner function 
thus does not become negative. However we have not proven that the initial classical 
state does not turn into a nonclassical state after photon addition and subtraction. 
Even though the negativity in the Wigner function is a strong proof of a state being 
nonclassical, the negativity or non-existence of the P function is another criterion for 
nonclassicality. 

In order to find out the nonclassicality of a state after adding and subtracting a 
photon for an initial classical state, we consider an example of the thermal field. After 
a straightforward calculation, we find the P function for p s B in ( 1921) : 



P s B a (a)=/V B a (n B + l)\a\ 



n B + l 2 
\ct\ 



e -|ar/n S) ( Q6 ) 



n B 

which becomes negative when the term in bracket is negative. Thus the resultant state 
is nonclassical. As a comparison, the P function for p B s in (19"2]) is found to be 



p a n s ( a )=m 



ran 
B V*) — ■ /V B 



riB + l, |2 ' 

\a\ — 1 



2 



n B + l: ,9 _|„|2 



n B 



a|He- |a ™ (97) 



n B 

The P functions are plotted in Fig. [7] for n B = 0.5. 

The quasiprobability functions have been considered for two initial states. Any 
sequence of photon addition and subtraction brings the state into the nonclassical state 
in our studies. In fact, the sequence of photon addition and subtraction has been 
experimentally realised for the input thermal field [31]. In this work, the reconstructed 
Wigner function shows the negative dip for the photon subtracted and added state. 

5. photon subtraction to increase entanglement 

5.1. A brief review of nonlocality, entanglement and quantum teleportation 

When two systems are correlated, knowing one system reduces the uncertainty of the 
other system. The correlation allowed by quantum mechanics is different from the 
conventional concept of correlation. The paradoxical concept of quantum correlation 
started from Einstein, Podolsky and Rosen (EPR) who questioned the completeness of 
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Figure 7. The P functions for p B s which is obtained after subtracting then adding 
a photon to a thermal field (a) and for p s g which is obtained after adding then 
subtracting a photon for a thermal field (b). The average photon number for the 
initial thermal field is ub = 0.5. 



quantum mechanics based on local realism in their seminal paper [56]. It seems that in 
this question, Einstein was not fully content with the indeterministic nature of quantum 
theory. A class of theories called local hidden-variable theories emerged in the search 
of a complete theory. About 30 years after the EPR paper, Bell [57] came up with an 
experimentally testable inequality which a local hidden-variable theory should satisfy. 
Clauser, Home, Shimony and Holt (CHSH) [58] 's version of Bell's inequality for a spin-| 
system is 

(AB) + (A'B) + (AB') - (A'B') < 2 (98) 

where A, A' are two observables, taking values either 1 or — 1, for one system and B, B' 
are two other observables for the other system. However, for the singlet state 



\tysingl 



^=(1 T>.| i>6 



l)a| ])l 



(99) 



the correlation function on the left hand side of the inequality can be as large as 2\/2, 
which violates the local hidden variable theory. This was experimentally proven using 
the two polarisation states of photons, instead of ups and downs of spins (See a recent 
experimental result [59] for some hints and problems in the experimental proof of the 
violation) . 

Differently from the spin system, which can have two measurement outcomes of 
spin up and spin down, the light field we consider in the paper is defined in an infinite 
dimensional Hilbert space. For instance, the number of photons in the field can be 
0, 1, 2, • • • , oo. In phase space, the light field is described by a continuous variable. 
There have been efforts to consider the nonlocality of light fields, using dichotomic [60] 
and multichotomic [61] observables as well as continuous- variable correlations [62] . 

Apart from the test of nonlocality based on the EPR paradox and Bell's inequality, 
the concept of entanglement has been extensively studied as one of the key resources 
for the current development of quantum information processing. While the test of 
nonlocality is using physical observables, entanglement started as a mathematical 
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concept. When a composite system of two subsystems a and b is not entangled 
(separable), the density p of the composite system can be written as a statistical sum 

p = J2p(i)p?Pi (100) 

i 

where p(i) is a probability function and p" and p\ are density operators for subsystems 
a and b, respectively. A density operator is a positive normalised operator and its 
transposition is again a density operator. Thus a partial transposition, 

/^ = X>(0l#] T & (ioi) 

i 

of a separable density operator ( 11001) has to bepositive, where the superscripts pt and 
t denote partial transposition and transposition^]- Thus, using the converse, we can say 
that the negativity of the partially transposed density operator is a sufficient condition 
for the entanglement of the composite system. The Horodecki family proved that the 
negative partial transposition (NPT) condition is also a necessary condition for the 
entanglement of certain classes of composite systems [63] • A few years later, Lee et al. 
proved that the NPT is an entanglement monotone thus can be used as a measure of 
entanglement for a spin-| system [51], which was extended to a high-dimensional system 
by Vidal and Werner [65] . For a Gaussian two-mode field, the NPT is also a sufficient 
and necessary condition [66] for entanglement. On the other hand, for a non-Gaussian 
two-mode field, a more elaborative entanglement condition applies [67]. 

Quantum entanglement is a mathematically well-defined concept and quite 
universal in the sense that entanglement in a composite system means its possession 
of quantum-mechanical correlation, regardless of its usefulness for a particular task. 
Quantum teleportation, on the other hand, is a physical task for which entanglement 
of a quantum channel is necessary for its success. Even though quantum entanglement 
is not a sufficient condition, the success of quantum teleportation can be regarded as 
a physically measurable witness of entanglement. In this sense, the nonlocality test 
discussed earlier is another experimentally accessible entanglement witness. 

After the two particles of a perfect entangled pair are shared between two remote 
stations, the quantum state to be teleported and one particle of the entangled pair are 
jointly measured at one station. The other particle of the entangled pair at the other 
station is then collapsed into a state. The original quantum information can be recovered 
by a unitary transformation according to the classical message on the measurement 
outcome (See Ref. [2] for a review on quantum teleportation of a coherent state and 
continuous variables.). 

+ A positive operator is defined as an operator whose eigenvalues are all positive. 

* For the purpose of an entanglement test or an entanglement measure, the partial transposition of 

either of the two modes will bear the same result. 
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5.2. Teleportation with a Photon- subtracted squeezed state 

For a two-mode (modes a and b) continuous-variable system defined in an infinite- 
dimensional Hilbert space, a maximally entangled state is 

oo 

\^max-un) = ^ \ n )a\ n )b (102) 

n=0 

which is unnormalised and the energy is infinite hence the state is unphysical. On 
the other hand, the two-mode squeezed state l^) in Eq. (TT9l is entangled with a finite 
degree. The degree of entanglement for a pure state can be obtained by the von Neumann 
entropy of the marginal density operator of a subsystem. For the two-mode squeezed 
state of density operator p\s 2 ), the degree of entanglement (£ is 

£(/5|s 2 >) = -Tr a [p a log 2 p a ] (103) 

where 

1 oo 

p a = Tr b (p lS2) ) = — ^tanh 2n CNH. (104) 

' n=0 

The marginal density operator of the two-mode squeezed state is already diagonalised 
so that the calculation is rather straightforward but for a general two-mode state, 
diagonalisation may be a challenging task. On the other hand, 97l(p) = 1 — Tr a p 2 , 
based on a linearised entropy, can be more easily calculated even though it does not 
provide statistical information as well as the von Neumann entropy does. 

Because of the experimental limitations, a highly entangled state is difficult to 
obtain. For example, in the first demonstration of the continuous-variable quantum 
teleportation [68], the squeezing was low hence the low entanglement of the quantum 
channel and the average fidelity of quantum teleportation was 58 ± 2%, which is 
just above the classical limit of 50%. In order to improve the fidelity of quantum 
teleportation, Opatrny, Kurizki and Welsch (69] suggested to select an optimal 
subensemble of entangled fields based on a conditional measurement, which is in effect 
subtracting a photon from each mode of the two-mode squeezed state. Once the 
quantum channel is prepared in this way, the rest of the quantum teleportation protocol 
is performed. For the initial state of a coherent superposition state, they show the 
improvement of the teleportation fidelity. At the end of their investigation, Opatrny and 
colleagues briefly discuss that they do not find an improvement of quantum teleportation 
for the quantum channel prepared by adding a definite number of photons to each 
mode of the two-mode squeezed state. Olivares, Paris and Bonifacio [70] . on the other 
hand, considered the photon subtraction using on-off photodetectors and showed the 
improvement of quantum teleportation depending on various parameters involved. 

The von Neumann entropy measure (£ of entanglement for the two-mode squeezed 
state when the squeezing parameter ( = 1 is 2.34 while that is increased to 3.53 
for the photon-subtracted squeezed state. The linearised entropy measure DJl(p) of 
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entanglement also shows that the photon-subtracted squeezed state is more highly 
entangled than the original squeezed state for non-zero squeezing: 

Tr a (Tr 6 p| 52> ) 2 >Tr a (pl s ) 2 (105) 

where p* s = Tibp ts and p ts is the density operator for the field after a photon is subtracted 
from each mode of the two-mode squeezed vacuum, 

1 oo 

P l : = cosh 4 CcQsh2C £(tanhC) 2 > + l)>>.<n|. (106) 

' ' n=0 

Cochrane, Ralph and Milburn [TTj calculated the photon number distribution P(n) for 
the marginal state p^ s : 

P{n) = (n\ft\n) = / tanh 2 " ((n + l) 2 (107) 

cosh (coshi(| 

As can be seen clearly in Fig. El the photon number distribution is peaked sharply at 

n = for the two-mode squeezed state while for the single-photon subtracted squeezed 

state, the peak is moved to n = 3 and the width of the distribution is wider with a 

smooth peak. 

Entanglement has two important ingredients: random and deterministic nature. 

As can be seen in the singlet state ( j99l) . before the measurement a subsystem is equally 

likely be measured in spin up or down; hence the random nature. However, once the 

measurement is performed on a subsystem, the state of the other subsystem collapses 

into an eigenstate of the measurement so that we can predict the measurement outcome 

without measurement; hence the deterministic nature. In order to further the analysis, 

let us consider a pure entangled state \ip ge n) — Cq\ T)o| l)b~ c i\ l)a\ T)& with real c and 

C\. When cq = c\, the outcome of the first measurement is maximally random. Now, 

we perform local unitary transformations U to rotate spins in such a way, 

IT>^^(IT> + li» , II)^-^(IT)-ID). 



W 9 en) -- ol !>4(CQ - Cl)| |) " (O) + Ci)| |)] 6 + -| |) a [(c + C X )\ ]) - (c - d)| |)] 6 . 



The \ipgen) turns into 

i|t>a[(cd-Ci)|t>-(cb + Ci)|i>]6 + i| 

Two states (c =F Ci)| t) — ( c o =t c i)| I) f° r particle b are not orthogonal to each other 
unless Cq = c 2 . Because of the ambiguity due to non-orthogonality, the correlation 
between particles a and b is not as clear as for the singlet case (\i^ gen ) when c = c\). 
Extending this idea, we can say that the enhancement of entanglement is related to the 
flatter weights of each component states in the photon subtracted squeezed statq$|. 

As discussed earlier in Eq. ( 1551) . the experimental realisation of single-photon 
measurements is approximated by an on-off photodetector. In this case, once a 
measurement is done on one mode, the field of the other mode is collapsed into a mixed 
state. Kitagawa et al. studied the degree of entanglement for the photon subtracted 
two-mode squeezed state with this experimental consideration [72] . 

t) This should not be generalised too much but at least this gives an intuitive view into the problem. 
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Figure 8. Photon number distributions for mode a of a two- mode squeezed state 
before (a) and after (b) subtracting a photon from each mode. Squeezing parameter 



5.3. Nonlocality of a photon- subtracted squeezed state 

While the original EPR paradox was on a continuous-variable system, Bell, CHSH and 
many followers developed inequalities for the so-called Bohm's version of the EPR state, 
which is for a spin-^ system. This system has only two measurement outcomes so that 
it is much easier to treat than the infinite-dimensional continuous variables. (There 
have been discussions on the validity of a nonlocality test for a multichotomic, many 
observables and many-dimensional system. As this is beyond the scope of the present 
paper, we will refer to Ref. [61] and references there in for further studies.) 

Bell wrote that the EPR state, which is simulated by a two-mode squeezed state 
when the squeezing parameter is infinity, would not violate his inequality because 
its Wigner function is positive everywhere in phase space. Using the fact that a 
Gaussian state transforms into a non-Gaussian state by subtracting a photon, Nha 
and Carmichael [73] suggested the CHSH Bell's inequality test for a two-mode squeezed 
state, conditioned on subtraction of a photon from each mode of the state. Their 
dichotomic observable is based on a homodyne measurement, which has high detection 
efficiency to avoid a detection loophole in an experimental nonlocality test. As shown in 
Eq. (l27j) . the phase value is measured by a homodyne detector. The observable setting 
is determined by the phase of a local oscillator (0^ in Eq. ( 1271) ). 

For the dichotomic value, +1 is endorsed when the measurement outcome is positive 
and -1 otherwise. In their analysis, Nha and Carmichael considered the realistic way 
to subtract a photon using a beam splitter and an on-off photodetector. They found 
that the CHSH Bell's inequality is violated even when the squeezing parameter ( is 
as small as 0.48 as far as the transmittivity of the beam splitter is high (~ 1). They 
also found that the scheme is quite insensitive to the detection efficiency of the on-off 
detector. It is because if a photon is not detected, they will anyway throw away the data 
so losing a photon is not a problem. In the same sense, losing a photon does not cause a 
problem in the on-off photodetector. Nha and Carmichael carefully remark against the 
possible fallacy due to the selection of data as the test is conditioned on subtraction of 
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a photon. They say [73] 'the situation differs from the so-called detection loophole which 
asserts that unfair data sampling can lead to the violation of a Bell inequality even for 
a classcially correlated state.' 

Garcia-Patron et al. [73] also studied the same system and assert that the local 
hidden variable theory does not change at all by the subtraction process because it 
cannot be influenced by the two local homodyne measurements. They considered the 
detection efficiency of the homodyne detectors and found 1% of CHSH Bell's inequality 
violation for the homodyne efficiency of about 95%. In a longer version of this work, 
Garcia-Patron, Fiurasek and Cerf [71] extended the scheme to consider various scenarios 
to subtract photons using beam splitters and on-off photodetectors. They found that 
the violation would increase by subtracting one, two and four photons from each mode 
but no violation occurs by subtracting three photons. 

5.4- Experiment: entanglement enhancement by single counting 

Despite theoretical studies on increasing nonlocality violation and teleportation fidelity, 
experimental realisation has not been made to subtract photons from each mode of the 
field. It is because the probability of subtracting a photon from one mode is already 
very rare and difficult so that coincidence subtraction of two photons is technically 
extremely demanding. However, there is a scheme which is less demanding because 
only one photon subtraction is involved to increase entanglement. If the subtraction is 
ideal for the operation to be represented simply by an annihilation operator, a single- 
photon subtraction from one mode of the two-mode squeezed state will increase its 
entanglement. The two-mode squeezed state turns into 

2 °° 
\S 2s ) = —^z J2 Vn(- tanhC)> - l)«|n)&. (108) 

n=l 

The linearised entropy entanglement measure for this is 

M(p 2s ) = 1 - (-r-^ J £ n 2 tanh 4 - C (109) 

which is higher than its counter part for the initial two-mode squeezed state for all 
nonzero squeezing parameter (. 

Ourjoumtsev et al. [75] use a technique to subtract only one photon but by erasing 
the information on the path information of the subtracted photon, they increase the 
degree of entanglement. Let us consider a scheme to annihilate a photon from one of 
the modes of a two-mode field. Then the operation is described by (a + b)/y/2. After 
this operation, the two-mode squeezed state \S 2 ) turns into 

/2 °° — 
\S 2e ) = -^— J2 ^(- tanhC) n (|n) a |n - 1) 6 + \n- 1)» 6 ) (110) 

The marginal density operator for mode a can be calculated as p a = Tr6|<S 2e )(<S2 e |: 



p a = ^ — > ntanh 2 " C 

H sinh 2 2( ^ 



n=\ 



\ n )\ n \ + i — r^\/ \ n + l/v 1 - i 

tanhC V n 
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Figure 9. Ourjoumtsev et al. 
[75] increases the degree of entanglement by subtracting one photon from one of the 
mode of a two-mode squeezed state. PD: photodetector, and T: transmittance of the 

beam splitter. 
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111) 



whose OJl is larger than VJl(p\s 2 ))- 

With the experimental setup as in Fig. [9], Ourjoumtsev et al. show an enhancement 
of entanglement for small squeezing. In the experiment, as an on-off photodetector is 
used, the resultant state in the experiment is mixed. Thus the von Neumann entropy 
or the linear entropy cannot be used to measure the entanglement. They reconstruct 
the density matrix from the homodyne measurements, then using the entanglement 
monotone based on the NPT condition [65] they show that the photon-subtracted state 
has higher entanglement 



6. Final remarks 

Since its advent, a laser light field has been an important tool to provide experimental 
evidences of paradoxical ideas in quantum physics. In the current development of 
quantum information processing and quantum control of a system, it is very useful to 
generate a light field in an arbitrary quantum state. In this paper, we have revised how 
to engineer the quantum state of a travelling light field using beam splitters, parametric 
amplifiers, photodetectors and homodyne detectors. We have seen that single-photon 
subtraction and addition are possible and these have been demonstrated experimentally. 
It is also possible to increase entanglement of a two-mode light field using the same 
technique. 

Even though a sequence of subtraction or addition enables to engineer a quantum 
state to an arbitrary state, it is experimentally demanding. It will be interesting to 
have theoretical suggestions which may be easily realisable by experiments for tests of 
quantum physics and to increase the applicability of quantum optics. 
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